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Abstract 

We derive a closed form portfolio optimization rule for an investor who is diffident about 
mean return and volatility estimates, and has a CRRA utility. The novelty is that confidence 
is here represented using ellipsoidal uncertainty sets for the drift, given a volatility realization. 
This specification affords a simple and concise analysis, as the optimal portfolio allocation 
policy is shaped by a rescaled market Sharpe ratio, computed under the worst case volatility. 
The result is based on a max-min Hamilton-Jacobi-Bellman-Isaacs PDE, which extends the 
classical Merton problem and reverts to it for an ambiguity-neutral investor. 
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1 Introduction 

Traditionally, financial modelling heavily relies on the choice of an underlying probability mea¬ 
sure P, which is chosen to incorporate the statistical and stochastic nature of market price 
movements. As early back as the works of Bachelier, Samuelson and Black, Scholes and Merton, 
the underlying risk factors—such as stock prices or interest rates—have been modeled as Marko¬ 
vian diffusions (with possible jumps) under P. However, as has become quite agreed upon, the 
complexity of the global economic and financial dynamics render impossible the precise identifi¬ 
cation of the probability law of the evolution of the risk factors. Unavoidably, financial modelling 
is inherently subject to model uncertainty, which also appears under the appellation of Knight¬ 
ian uncertainty. 
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mustafap@bilkent.edu.tr. Department of Industrial Engineering Bilkent University 06800 Bilkent, Ankara, Turkey. 
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In the presence of model uncertainty, one may admit various degrees of severity. One may 
deal with model misspecification only at the level of the equivalence class of P, or go beyond and 
take into account a family of non dominated models. The core issue of portfolio optimization has 
been widely investigated over the last twenty years in the multiple priors context. The investor 
has a pessimistic view of the odds, and takes a max-min (also known as robust) approach to the 
problem, first minimizing a utility functional over the priors and afterwards maximizing over the 
investment strategies. We are aware of only a few results in the non dominated case, notably 
Hernandez-Hernandez and Schied [9] and the recent preprints by Nutz m in full generality but 
discrete time, and Lin and Riedel m in a diffusion context. On the contrary, in the dominated 
priors case there is a rich literature. We content ourselves with citing Chen and Epstein [3], 
Garlappi et alii HH, Maenhout Follmer et alii HH! for a comprehensive review and refer¬ 
ences, and, more recently, the work by Owari m- 

In such an active environment, the present note offers a resolution of the robust non- 
dominated Merton problem, which is both simple and mathematically rigorous. The main 
novelties of the present contribution lie in the form of the uncertainty set and in the accommo¬ 
dation for market incompleteness. We assume that the asset prices process is an N -dimensional 
diffusion, and the driving Wiener process is d-dimensional with d > N. The investor is diffident 
about the constant drift and volatility estimates /I and a. Thus, she considers as plausible all 
the variance-covariance matrices lying in a given compact set and, for a given realization of a, 
she considers all the drifts which take values in a ellipsoid centered at jl : 

Ue{a) = {u G I (u — — ft) < e^}, 

in which e > 0 is the radius of ambiguity and S = aa' is the variance-covariance matrix. 

The merits of an ellipsoidal representation for the ambiguous drifts has been amply demon¬ 
strated and discussed in m, m for the robust mean-variance optimization. The problem of 
worst-case (max-min) robust portfolio choice is a well-studied problem (see e.g., H El El El El 
EH Ea EZl EB] for robust portfolio optimization in single period problems) under different rep¬ 
resentations of ambiguity. Intuitively, the non-linear but simple geometry of ellipsoids offers 
robustness that avoids a worst case which is a corner solution. This is the case in a polyhedral 
hyper-rectangle or box representation, as in Lin and Riedel where the drift (as well as volatility) 
is allowed to vary in a box In the dominated setup, the assumption of A:-ignorance in [3] 

also amounts to a box representation for the drift. 

At the same time, the choice of U^^cr) preserves tractability. Citing Fabozzi et al. [6]: “The 
coefficient realizations are assumed to be close to the forecasts, but they may deviate. They 
are more likely to deviate from their (instantaneous) means if their variability (measured by 
their standard deviation) is higher, so deviations from the mean are scaled by the inverse of 
the covariance matrix of the uncertain coefficients. The parameter e corresponds to the overall 
amount of scaled deviations of the realized returns from the forecasts against which the investor 
would like to be protected.” 
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Another appealing feature of taking model uncertainty into account is that it offers a theo¬ 
retical solution to the equity premium puzzle. As noted by Mehra and Prescott [Ej, the high 
levels of historical equity premium and the simultaneous moderate equity demand seem to be 
implied by unreasonable levels of risk aversion. Their conclusion was skeptical on the ability of a 
frictionless Arrow-Debreu economy to account for such empirical evidence. However, the works 
by Abel [1] and Cecchetti, Lam and Mark [2] addressed the equity premium puzzle by relaxing 
the hypothesis that the investor perfectly knows the probability law. The key point is that, in 
the multiple priors setup, the optimal equity demand depends on two aversion components: risk 
and ambiguity aversion. In accordance to these results, and the subsequent [3], [13] and m, 
we find that robustness of decisions lowers the optimal demand of equity since the ambiguity 
and risk averse investor effectively behaves like a risk averse investor with an increased risk 
aversion coefficient. Precisely, in a CRRA utility case with relative risk aversion parameter R, 
the optimal relative portfolio is given by 




iH-e)+ 

RH 




— rl) 


in which S is the worst case variance-covariance matrix, and H is the Sharpe ratio computed 
under S (see Proposition 1, 2 and Section 4). So, when the ambiguity radius is too high, namely 
it exceeds the worst case Sharpe ratio, the investor refrains from investing in the risky assets 
and puts all the money in the safe asset. The opposite case is when there is no uncertainty 
in the drift, that is e = 0, and no uncertainty on volatility as well: then the optimal solution 
reverts to the Merton relative portfolio, ttm ■= ~ ^1)- 

The paper is organized as follows. Section 2 contains the (diffusion) model specifications in 
the non dominated case, together with the general version of the Martingale Principle needed 
here. However, to derive an abstract max-min PDE from the Martingale Principle, some con¬ 
ditions on the volatility structure must be imposed in the fully incomplete market model, as 
in |9] for the case of a traded asset with coefficients depending on an underlying, non traded, 
asset. The focus here is on the complete market case, to wit the volatility is a square matrix. 
This allows for a simpler, yet effective, analysis. In fact, the HJB-Isaacs PDE formulation shows 
that the investor is observationally equivalent to one who has distorted, worst case, beliefs on 
the parameters. In Section 3, a representative investor with CCRA utility is considered, and we 
assume that ambiguity is present only in the drift, i.e. the priors are all equivalent. This is an 
interesting case per se, since the drift is subject to imprecision in estimations to a much greater 
extent than volatility. There, we solve and provide the explicit solutions to the robust problems 
for the infinite and finite horizon planning. Finally, we apply these findings in Section 4 to some 
examples in the genuinely non dominated setup. 
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2 The Merton problem under ambiguity aversion 


Consider the problem of an agent investing in n risky assets and a riskless asset. Specifically, 
we work under the Black-Scholes-Merton market model assumptions. Namely, the riskless rate 
r is constant and the risky assets dynamics are, for each i = 1,... ,d: 


dSi = si f/dt + 


d 


( 1 ) 


where and are constants and W is a standard, d-dimensional Brownian motion on a 
filtered space (0, (iFt)t>o, P) ■ Assume that d > n, so the market is allowed to be incomplete. 
In matrix-vector form, the above equation becomes: 


dSt = Diag{St){fJ.dt + adWt) 

where by Diag{St) we denote the diagonal n x n matrix with i-th diagonal element equal to 
SI. In addition, a is required to have full rank, so that the variance-covariance matrix S = aa' 
is invertible. Here and in what follows, ' denotes the transpose operation. 

Given the initial endowment wq , the investor is allowed to trade and consume in a self¬ 
financing way. To be explicit, let h = {ht)t denote the n-dimensional progressively measurable 
process, representing the number of shares of each asset held in portfolio, and let the progres¬ 
sively measurable, nonnegative, scalar process c indicate the consumption stream. Assume also 
that Jq Cgds is finite P—a.s. Then, the wealth process is governed by the following stochastic 
differential equation: 

dwt = {rwt -|- hiDiag{St){g — rl) — ct)dt + h[Diag{St)(jdWt 


in which 1 is the d-vector with all components equal to one. It is convenient to recast the wealth 
equation by the vector process 6 of cash value allocated in each risky asset, i.e. 9t := Diag{St)ht ■ 
Thus, 

dwt = [rwt + Oi{p — rl) — ct)dt + O'adWt- (2) 

The pair {9t,ct) is admissible for the initial wealth wq if the wealth process wt given by (§ 
remains P-a.s. non-negative at all times. Let A^{wo) be the set of all admissible (d,c) pairs 
for initial wealth wq. Note that the admissible set depends only on the equivalence class of P. 
The agent is then trying to choose (d, c) G A^{wo), so as to maximize the expected utility from 
running consumption and terminal wealth: 

sup K[ f u{t,ct)dt + u{T,wt)]- 
{e,c)eAPiwo) Jo 

The utility function u is assumed to be concave, increasing in the second argument and measur¬ 
able in the first. This class of stochastic control problems is known under the name of Merton 
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problem. It includes a number of specific cases, among which the infinite horizon planning. 
In fact, if T = oo, Woo ■= lim sup^^^^ rcj and n(oo,-) = 0 the above optimization problem 
becomes: 

poo 

sup E[ / u{t,ct)dt]. 

{0,c)&A^ (wo) Jo 

So far, the exposition is classical, and can be found in many textbooks. The reader is referred 
to the new m Chapter 1], for a remarkably didactic approach. 

However, things change quite a bit if the agent is diffident about the (constant) estimates 
jl and full rank matrix d, for the drift and volatility matrix of the risky assets respectively. 
Assume from now on that II is the Wiener space of continuous functions, with the natural 
filtration F = Our investor assumes that the ‘true’ volatility a is a progressively 

measurable matrix, and such that the variance-covariance S = aa' takes values in some fixed 
compact set K of n x n invertible matrices, containing S: 

S-.= {a e I S e AT}. 

Let us denote by 5 = {cr progr mis | at{oj) G S for all ca, t}. This choice is in line with empirical 
practice, as S is the estimated object, not the volatility a. The uncertain drift is also assumed 
to be progressively measurable, and for a given realization of a it is allowed to vary in 

Ue{a) = {u G M"' I (u — — fi) < for all w,!}, 


that is, in an ellipsoid centered at jl with radius e. Let us denote the set of plausible drifts and 
volatilities by 

T := {(//, fj) progr. meas. | a G G Ue{at{u:))] 


and let denote its ci-section. Different choices of (/i,cr) G T correspond to considering 
different probabilities on the Wiener space, namely those under which the risky assets evolve 
with the prescribed coefficients. These probabilities are orthogonal to each other across the 
sections To-. 

For a fixed choice of the process a however, the Girsanov theorem ensures that all the vector 
processes fi G correspond to probabilities that are equivalent to each other on T) for all 
t > 0. We would like to describe these equivalent changes of measure as a function of /U. To 
this end, let us select a reference probability corresponding to ( jlyCr) and call it . Market 
incompleteness implies that the probability in under which the risky assets evolve with 
drift /r is not unique. However, such measures can be fully parametrized as probability changes 
with respect to . And a minimal choice (see also Remark below) is selecting for each 
/i the probability P^^’^ corresponding to the measure change given by the Doleans exponential 
T(J’q (ff^dW ), where 


Tt' := Hd-t - jl)- 
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Such selection does not reduce generality, since what matters in our context of expected utility 
maximization are only the distributional properties of the risky assets. Therefore, we have now 
a one-to-one correspondence between elements of T and probabilities on (n,F), namely the 
possible prior models are given by 


Now, the wealth process evolves under each according to 

dwt = (rwt + 6[{nt - rl) - ct)dt + 6[atdW 

where IT is a -standard Brownian motion. Finally, let us call the investment/consumption 
pair (9, c) (robust) admissible for the initial positive wealth wq if in addition to the measurability 
and integrability assumptions already made at the beginning of this section, the wealth process 
remains non-negative for all P G P; 

>^rob('fi’o) := (wo) = (w^o)- 


The equality on the rhs holds because, for a given a , the admissible class is invariant for different 
choices of /r G . The ambiguity averse investor takes a prudential worst case approach, and 
faces the following robust Merton problem: 


Uoptiwo) ■= sup inf 
(e,c)eAab(«'o) 


Ljo 


u{t, ct)dt + u{T, wt) 


( 3 ) 


It is clear that more conservative portfolio choices are made when the uncertainty set T is larger, 
while an ambiguity-neutral investor sets S = {d} and e equal to zero, thus facing a classical 
Merton problem for the model . 


Remark 1 Fix a £ S. Any change of drift, say from fi to fi, corresponds to a change of 
measure from Ph’^ to a probability P with density process ^ given by a Doleans exponential 
i£’(Jq (fdW). The suitable ips can be characterized as those in the form: 

= + ijjt 

in which ip is a (sufficiently integrable) progressively measurable process, belonging dt® dP^'^ - 
a.e. to ker((Ti(a;)). Chen and Epstein call the process (p the market price of ambiguity. 
Denote this class of probabilities by . Elementary optimization shows (p^ is minimal, in the 
sense that it has the smallest pointwise M"" -norm among : 

||(/?^(ti;)||^ = min ||(^t(a;)|p dt ® dP^’^ — a.e. 

Pg-P<T 

Thus, the ellipsoidal ambiguity on the drift in can be recast into a (non convex!) condition 
on the market price of uncertainty (p, namely 

/r G To- iff min \\(pt{(jj)\(^ < e^ dt ® dP^'^ — a.e. 

Pg-P<T 
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This should be contrasted with the ubiquitous conditions in the literature (see e.g., m,m) which 
require the market price of uncertainty to be valued in a convex set. 

The resolution of Q is based on the next robust version of the verification theorem. Although 
formulated for the sets V,T, its validity is general and does not rely on any specific parametriza- 
tion of the set of prior models. 

Theorem 1 For a shorthand, call v = {p,, a) a generic element of T. Suppose that: 

1. there exists a function V : [0,T] x —>■ M, which is continuous on [0,T] x M"*“ and 

on [0,T) X M+, verifying V{T, .) = u{T, .); 

2. for any [Q, c) there exists an optimal solution n{6, c) G T of the inner minimization in 
([sl), such that 

Yt = = V(t, wt) + / u{s, Cs)ds (4) 

Jo 

is a -supermartingale; 

3. there exist some {0,c) G ^rob('«^o) such that the corresponding Y is a - martingale. 

Then {9, c) is optimal for the problem (© and 14(0, rro) is the optimal value function, namely 

Uopt{wo) = 14(0, rco) • 

Proof: The proof is a simple modification of the Davis-Varaiya Martingale Principle of Optimal 
Control |20l Theorem 1.1]. In fact, by the supermartingale property of Y under ^ and by 

14(r,.) = u{T, .), we have: 

^ ^U{e,c) J ^ ^ 

40 

Taking the supremum over Aob(ioo) gives Uopt{wo) = sup_ 4 ^^j^(^p) n(s, Cs)ds-|-u(T, ict)] < 

14(0, Wo). Since by assumption for some V = {6, c) the process T is a martingale under , 

then K'^[Yj']=Yo = 14(0, rco) and the conclusions immediately follow. | 


Now, the usage of the verification theorem to solve the ambiguity-averse investor’s problem 
is quite intuitive. Given a specific utility function, one looks for a function 14 satisfying the 
premises of the theorem. Using ltd’s formula, any process Y as in (Q verifies under the 

following SDE: 

dYt = {u{t, ct) + Vt-\- Vw{rwt + 9[{iJ.t - rl) - c*) -F ^etYt6tV.ujw}dt + VwO'tcrtdW. (5) 

To make Y a supermartingale under every ^ and a martingale for some ^ the 

maximum over {0, c) G M” x of the minimum of G T must be equated to zero. At this 
point, some other specific structure on a must be assumed in the fully incomplete market case, 
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like e.g., dependence on a correlated nontraded asset as in [9]. 

In the rest of the paper however, we focus on the complete market case, i.e. we assume a is 
a square matrix. Then, a max-min nonlinear PDE arises from ([^: 


max mm 

(0,c)£lR”'XR+ CrGS',/.iGf/e((T) 


u{t, c) + Vt + Vw{rw + 9'{n- rl) - c) + -9'T,9V^ 


= 0 , 


which is of the Hamilton-Jacobi-Bellman-Isaacs type. In the following, we simply refer to it as 
the robust HJB equation. The minimization of the drifts for a fixed gives: 


mm.{9'u : (u — llYT, ^(u — n) < e^}. 


which is a simple exercise in Karush-Kuhn-Tucker necessary and sufficient conditions. 
0 7 ^ 0, the optimal solution is 


^{9) =fi-e 


VWYY9' 


When 


Substituting it back in the robust HJB, we get 


max mm 

(6,c)sIR"xR+ crS5 


U 


{t, c) + Vt + Vw{rw + 9'{fi - rl) - eV9'T.9 - c) -|- ^9'T.9Vw 


= 0 , ( 6 ) 


which covers also the case 0 = 0. Since the value function V will be increasing (and concave) 
in the wealth w , we are minimizing a concave function over S: 


min 

ctGS 






(7) 


Notice that the function to be optimized depends on a only via the quadratic form 9'T,9, and 
that its derivative wrt y := 9'T,9 is positive as 14; > 0, < 0. Therefore, the optimizers are 

not unique in general and compactness of S is of essence. The actual computations depend on 
the specification of 5, and may be quite easy when the set S is defined precisely in terms of 
constraints on the quadratic form, as we show in Section]^ 

Let us denote by d'{9) an optimizer, where d'(0) is the Cholesky factorization of an optimal 
varcov matrix S(0). Since the optimal value of quadratic form 9'T,{9)9 does not depend on the 
specific choice of d'(0), one obtains the PDE 


max 

(6»,c)eIR"xR+ 


u{t, c) + Vt + Vu,{rw -I- 9'{fi - rl) - e\J9'T^{9)9 - c) + ^0'S(0)0I4) 


= 0 , ( 8 ) 


In the main applications we present in Section 3 and 4, E(0) = S, namely a constant. When 
this is the case, the above equation is equivalently viewed as stemming from the worst-case 
(/i,fT), in which a is the Cholesky factorization of S, y = fj, — e-^=, and using the worst 
couple in the wealth equation: 


dwt = 


{rwt + 9[{fit - rl) - e\J9'YL9t 
8 


ct)dt + 9[adW. 


(9) 
















Therefore, the general problem Q becomes equivalent to a robust utility maximization when 
there is (ellipsoidal) uncertainty in the drift only. 

The techniques employed to solve Q are then standard, and rely on educated guesses at the 
form of the solution. If the solution V of the robust HJB can be found explicitly, then it is the 
candidate to be the value function we are looking for. Finally, to conclude that V is indeed the 
value function, one must check that it verifies items 1, 2 and 3 in Theorem 


3 The robust power utility problem with non ambiguous a 

We assume in this Section that the investor has a power utility function, and there that is no 
uncertainty on the (constant) square volatility matrix, namely S = {d} . Lack of uncertainty in 
the volatility may be empirically justified by the consideration that mean returns are subject to 
imprecision to a much higher extent than volatilities. We provide explicit solutions both when 
the planning horizon is finite and infinite. To avoid notation overload, and for the next usage 
in Section 4, we drop the hat over a. Also, we denote by 

lAe := {p progr meas | /rt(w) G C^e(o') for all w}, 
and by the expectation under . 


3.1 The infinite horizon planning 

3.1.1 Resolution of the robust HJB equation 


Let us assume the investor has CRRA power utility from intertemporal consumption: 

p-R 

u{t,x) = 

where p and 1 are positive constants, modeling the time impatience rate and relative risk 
aversion respectively. In the infinite horizon case, we wish to find the solution of: 


Uopt{wo,e) 


sup inf E'^ 



e 


-ps 


p-R 

'-s 

1-R 



( 10 ) 


when the problem is well-posed, i.e. when it has a finite valu^ Assume for the moment that 
this is the case and also that both the inner infimum (for a fixed {9,c) G .4.TOb('fi’o )) and the outer 
supremum are attained. The properties of the problem imply, exactly as in the classic case, that 
a guess at the value function takes the form 


V{t,w) = % %{t,w). 

^When S = {c}, we remark that Glrob(wo) coincides with the classic set of admissible plans (wq). 
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The positive constant 7 ^ has to be determined, and we use e as subscript to highlight the 
dependence on the radius of ambiguity e. With this guess, let us solve Q. The optimization 
over c trivially results in 

c = y.tc. 


with 

max{u(t, c) - cVw} = e~P \ ^ 
c \ — K 

The residual optimization is 


max 

6 


+ v^{rw + e'ift - ri) - + \e'T.ev^ 


l-R 


The function to be maximized is concave in 9, and smooth in \ {0}. The first order 
conditions are thus necessary and sufficient for optimality in 0 7 ^ 0. So, by equating the gradient 
to zero we obtain: 


e{s) = 


SVin 






where s := \^6'Ti6. We are left with 


= 0 (s)'S 0 (s) 


Set 

H := v^(/i-rl)'S-i(/l-rl), H, := H - e. 
The above equation has a positive root, given by; 




s = -- 


V,„ 


if and only if ifg > 0. If < 0, the optimal solution is necessarily 9 = 0. Finally, if 

denotes the positive part of , the following is a compact way of writing the optimal solution 

in both cases: 

9 = w^Y.-\fL-rl) 

Now, 7 e is found by substituting these c and 9 back into 
Straightforward calculations result in: 


and solving for the constant. 


7e = 


p+{R-l){r+li^) 


R 


( 11 ) 


1 h'^ 

which for e = 0 falls back to the constant 70 = ^ ^ classic case. Therefore, 

the value function V of the problem is found as 

V{t,w) = 


This of course holds as long as 7 ^ > 0, which is shown below to be a necessary and sufficient 
condition for the well posedness of the robust Merton problem. 
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3.1.2 The verification and comparison with the classic Merton problem 

Proposition 1 The infinite-horizon robust Merton problem under ellipsoidal ambiguity of mean 
returns: 


Uopt{'Wo,€)= sup inf 
(0,c)eAob(«>o)^e“' 


I 


oo „1-R 

.-psfis - 


l-R 


is well posed if and only if 7 ^ in © is strietly positive. In this case, the optimal value is 


Uopt{wo,e) = ^(OjWo) = 7e 


l-R 


and the optimal controls are: 

6t = WtTTe , ct = %Wt, 

with optimal portfolio proportions vector given by 


Hfi 


TT. := 


RH 




The worst case drift is constant: 


p.:= p{e) = fi-e 


\/ < Svr, 


:7r. 


and the optimal wealth process has dynamics given by 

( Tx. / 1) 

Wt = Wo exp ( TTeaWt + (r H-- -fpt 

Proof: The proof is split into two steps. 


( 12 ) 


1- If 7 e < 0 then Uopt{wo, e) = 00 . Note first that this case can only happen when 0 < i? < 1. 
The proof here closely follows the lines of eqi Section 1.6]. 

1-a) Assume 7 ^ < 0. Then, consider a couple of controls that are both proportional to 
the wealth: 

6t = wtT:, ct = Xwt with A > 0. 

If we substitute them into (7), then the solution is the positive wealth 


Wt = Wo exp 


^(r + 7r'(fi — rl) — eV vr'Evr — A — - 7 r'S 7 r)t + vrVlTt^ 


so that: 


Uopt{wo,e) > 



00 


/O 1 — i? 

An application of the stochastic Fubini’s theorem shows that the latter is proportional 
to 


exp 


t ( —p + (1 — ii)(r + 7 r'(/i — rl) — eVvr'STT — A — —tt'Ett) 


dt. 
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If there exist some (vr, A) for which the exponent is positive, then the integral diverges 
and the value function is infinite. For fixed A, the maximum over vr in the exponent 
is attained for vr = vTe, and the value is 

-p + (1 - R){r + - A) = -Rj, - A(1 - R), 

which is positive for A small enough. 

1-b) If 7 e = 0, take 9t = wtTTe, and c* = for some constant A; > 0. This choice 

leads to 

fOO 1 

Uopt{wo,€)> 

Now, when t —)■ 

oo. Therefore, the integrand is asymptotic to and hence the integral 

diverges if e.g., k = — 1. 


2. If 7 e > 0, then the optimal processes/value are as given in the statement of the proposition. 
Substituting the candidate optimal controls {9t, ct) into ([^, and solving for the candidate 
optimal wealth one gets (12), which can be further simplified to: 


wt = wq exp 


T^eWt + - p- 

K 



The process rc is a deterministic scaling of Geometric Brownian motion, as well as the 
process vt := Now, 


{1 - R)Yt := {1 - R)[V{t,wt) + f u(s,Cs)ds] = 7 g P^vt+ f e ^^( 7 e)^ ^Vsds. 

Jo Jo 

has integrable maximal functional in every compact [0, T] and by construction has zero 
drift term. Henceforth, T is a martingale. 

For what concerns other admissible controls {9t,ct), the process 


Yt = V{t,wt) + / u{s,Cs)ds 


with w as in (§, is by construction a diffusion, which has the same sign as (1 — i?), and 
which has non positive drift under . 


• If 0 < i? < 1, then any such Y is positive. By writing ([^ under it is imme¬ 

diate to realize that T is a positive, decreasing scaling of a positive local martingale, 
namely the Doleans exponential of the process X defined by: 


X := 


r p-ps 

/ ^:^—w-^9',adW. 
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Henceforth, H is a -supermartingale. In addition, I/(oo,-) = u(oo,-) = 0 so 

that the conditions of the Verification Theorem are satished, and the proof in this 
case is complete. 

• If P > 1, any V is negative. A simple modification of the argument just used in 
the 0 < P < 1 case only shows that V is a local supermartingale. Therefore in 
this case we show the optimality of {{6,c),fL) in another way. To this end, note that 
the martingale property of Y as above, together with standard minimax inequalities, 
gives 

E^[ J e~P^ “ Uopt{wo, e) < 

< inf sup EP[ e-P^ ^ dt] < inf sup E^[ / ^ dt], 

AlSWe _4^ob(uio) Jo 1 — P Arob(wo) Jo 1 — P 


SO if we prove that the first value on the left is equal to the value of the last problem 
on the RHS, we are done. This is quite an easy task. In fact, for a fixed constant 
^ G f/g the inner supremum is a standard Merton problem. Hence, 


sup E^ 

(6»,c)eTrob('i.«0) 



e-pt 


(Ci)^ ^ 

1-R 


dt 


(7(h))-'' 


(rco)^ ^ 
1-P 


in which we pose 7 (/i) := ^+1^ — ~ —-^(h) •= \/(h “ — rl). 

The residual minimization: 


int (bM)-" 


(wq)^ ^ 
1-P 


is then a simple exercise, the minimizer being //, so that 7 (h) = 7(h) = 7e) which 
ends the proof. 


I 


Let us remark that the optimal portfolio 0 preserves the form of the Merton’s Mutual Fund 
theorem. In fact, the optimal portfolio consists of an allocation between two fixed mutual funds, 
namely the riskless asset and the fund of risky assets given by — rl). At each time point 

the optimal relative allocation of wealth is now dependent on the ambiguity aversion of the 
investor in addition to his/her risk aversion through the coefficient: 

Hi 

RH' 

The above allocation naturally collapses to the Merton allocation ^ for e = 0. In case the radius 
of ambiguity e is greater than or equal to the market Sharpe ratio H, the optimal control policy 
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is not to invest at all into the risky assets. Since ^, the robust Merton portfolio TTg 

has smaller positions in absolute value with respect to the classical Merton portfolio. To wit, 
both long and short positions are shrunk with respect to the ambiguity-neutral portfolio. As 
expected, and already anticipated in the Introduction, robustness in the decisions lowers the 
optimal demand on equity, and thus offer a theoretical basis for a possible explanation of the 
equity premium puzzle. 

The consumption in the ambiguity averse case may be increased or curtailed, depending on 
the sign of i? — 1. In fact, when the classical problem and its ambiguity averse counterpart 
with e > 0 are both well posed, if 0 < ii < 1 then 7 ^ > 70 > 0 , while if ii > 1 the opposite 
inequality chain holds. 


3.2 The finite horizon planning for non ambiguous a 


Now the investor has a CRRA power utility both from intertemporal and terminal consumption 
at time T < 00 : 


u{t, w) = e 


= 


w 


l-R 


for 0 < t < T and u{T, w) = A 


w 


l-R 


l-R - ^ ' l-R 

in which A is a fixed positive constant. Here, we set the deterministic scaling of the CRRA 
power utility identical to that of the infinite horizon case to better highlight the similarities, but 
everything stated below holds also if is replaced by an integrable, positive and deterministic 
function h{t). We then wish to find the solution of: 


Uopt{wo,e)= sup inf 


r'T A-R 

^-ps±s - 


l-R 


ds A 


w, 


l-R 


l-R 


(13) 


Using the scaling properties of the CRRA utility, the guess to the value function is of the form 
V{t,w) = for some positive, differentiable function satisfying f{T) = A. The HJB 

equation (fsl) now looks like 


max 
(e,c)eiR"xE+ 




^{rw+9'{il—rl) — eV0'T,9 — c) — ^f{t)w ^ ^9'T,9 
1 — Ii 2 


l-R 

Proceeding exactly as in the previous section, one obtains 


= 0 . 


c{t, w) = w 


—pt \^/R 


(S) 


9 = WKf. 


Substituting the above back into the HJB equation results in a first order ODE for /: 

f f{t) + kj{t) + Re"fl*(/(t))^“s = 0 

1 f{T) = ^ 


with 


k, := (1 - R) -h <(/! - rl) - = {1 - R){r + 
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With the substitution f{t) = g{t)^, the ODE can be linearized and easily solved: 

g{t) = exp It R 

Comparing this to the solution of [201 Section 2.1] the only changes are: 1) the constant k^, 
in which is replaced by and 2) the optimal portfolio allocation, which is identical to 

the robust allocation case of the previous section. Obviously for an ambiguity neutral investor 
with e = 0 we fall back to the finite horizon solution of the Merton problem. 

Let us conclude by summing up the results just found, leaving the verification to the reader. 


Proposition 2 The finite horizon robust Merton problem under ellipsoidal ambiguity of mean 
returns 


Uopt{'Wo,€)= sup inf 
(0,c)eAob(«>o)^e^' 


e-ps^^ds + A^^ 


/o 


1-R 


l-R 


is always well-posed, and admits the optimal controls: 


= WtM 


ct = Wt 






where 


g{t) = Ar exp (^^{T - t)^ + e R^ exp ^ s^ds, 


and ke = {1 — R){r + ). The optimal fx = ia{9) = fi — e TTe, and the optimal wealth 

yTr'E-TTe 

process w has dynamics under given by: 


+ 'l2 , 


Wt = Wo exp 


t -R-s 
e R 




ds + tt' aWt 


4 Examples with ambiguous a 

In all the following examples, the volatilities are square, full rank, matrices. 

Example 1 (The uncorrelated case) Suppose that estimated volatility matrix a is diagonal. 
To wit, the risky assets returns are (instantaneously) uncorrelated. Further, we suppose that the 
ambiguity does not affect correlations, namely the ambiguity set S is that of diagonal matrices, 
whose diagonal S lies in some product 

[<d,al]x...[gi,al], 

with inf j CTj > 0 and < ai < ai. This is exactly the case examined by Lin and Riedel m, 
where the problem is treated via a G-Brownian motion technique. 
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For a fixed ambiguity radius e > 0 on the drift, the solution of the residual inner minimization 
over a in the max-min HJB becomes a triviality with this diagonal uncertainty specification. 
The unique worst case volatility is constant and it is the ‘highest’ one, d = Diagia\,... ,dn) 
and does not depend on 9. Therefore, the general problem (© becomes equivalent to a robust 
utility maximization with volatility a and ellipsoidal uncertainty on the drift only, with radius 
e. To give an explicit example, in the power utility case one ends up in solving (10) or (13) with 


a = a. It is clear then that the verifications are identical to the ones just seen in the previous 
Section. The resulting optimal relative portfolio is also constant: 


in which H = \J{ft — rl)'S ^{fi — rl),H^ = {H — e)+ 


Example 2 (Upper bound on the quadratic form S) This example can be seen as relax¬ 
ation of the previous one, in the sense that we do not impose constraints separately on each of 

the eigenvalues of T,, nor we assume that S is diagonal. We simply restrict the quadratic form 

—2 — 

induced by S not to exceed a given threshold A > 0 on the unit sphere, with A > \m , the 
latter being the maximum eigenvalue of d. This amounts to imposing the same bound on the 
maximum eigenvalue of S. Precisely, the volatility is assumed to be valued in 

S:={ae I 0 < x'S® < A^||xf for all x e / 0} 


The minimizers d of the inner minimization in ([^ are volatilities a with maximum eigenvalue 
equal to X, and such that 6 is an eigenvector relative to A. Finally, the robust HJB boils down 
to the concave maximization: 


max 

(6l,c)eR"xR+ 


u{t,c) + Vt + 14) (r 


= 0 


Therefore, the ambiguous volatility problem is observationally equivalent to a Merton problem 
with volatility matrix equal to XI, and drift uncertainty radius e. So, from here one proceeds 
exactly as in Section 3. The optimal relative portfolio is thus 


Trfia) 


Mu-rl\\ - 


fllA-^i|| A 


1 

-oiT-rl) 


(||/l-rl|| - Ae)+ 1 

^IIA-’’1|| A^ 


ifi-rl) 


Remark 2 Another interesting case of ambiguity specification on volatility is based on defining 
a ball around an estimate of the variance/covariance matrix. However, it does not lead to a 
closed form portfolio rule. 
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The ambiguity on volatility is represented here as membership to the set S' = {S ^ 0 : 
||S—S||i? < 5} for some estimate S of the variance/covariance and a positive parameter 

5. The inner optimization problem max^es V^e\/0'S0 — ^^9'Ti9 is equivalently posed as the 
matrix optimization problem in the space of symmetric n x n matrices: 

max V^e^/W^) " ^{99', S) 
ses 2 

where {X,Y) = Tr{X'^Y) is the trace product (inner product in the space of symmetric nxn 
matrices). Denoting rank-one matrix 99' (obtained from the dyadic product) as 0, and passing 
to variables X = T, — T, we have the equivalent quadratically constrained optimization (convex) 
problem in matrix variable X: 

max (0, X) + (0, S) + a{Q,X) 

subject to 

{X, X) < 6^ 

where we defined a = jgj- convenience and we omitted momentarily the 

term a(0, S). 

From the first-order optimality conditions, after some straightforward algebra we obtain the 
following optimal (worst-case) matrix 

E = E + i(a+|)e 

where A = (0, S) and B = ||0||fj f = VA + 6B, A = i ^^.at S is positive 

definite. Unfortunately, substitution of the worst-case matrix into the robust HJB equation 0 
yields a fourth-order polynomial which does not lead to a closed-form portfolio rule. Nonetheless, 
a numerical procedure can be used to find the optimal portfolio rule. 
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